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Continuing our investigation of the regularization of the noise kernel in curved spacetimes [N. G. 
Phillips and B. L. Hu, Phys. Rev. D 63, 104001 (2001)] we adopt the modified point separation 
scheme for the class of optical spacetimes using the Gaussian approximation for the Green functions 
- - - a la Bekenstein-Parlcer-Page. In the first example we derive the regularized noise kernel for a thermal 

field in fiat space. It is useful for black hole nucleation considerations. In the second example of an 
, optical Schwarzschild spacetime we obtain a finite expression for the noise kernel at the horizon and 

recover the hot flat space result at infinity. Knowledge of the noise kernel is essential for studying 
issues related to black hole horizon fluctuations and Hawking radiation backreaction. We show that 
^ |. the Gaussian approximated Green function which works surprisingly well for the stress tensor at the 

D ■ Schwarzschild horizon produces significant error in the noise kernel there. We identify the failure as 

^0 ' occurring at the fourth covariant derivative order. 

I , PACS numbers: 04.62.+V 
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ON ■ The noise kernel is the vacuum expectation value of the stress-energy bi-tensor for a quantum field. In curved 
. spacetime field theory [|l|, |[ |[ |4| it plays a role in stochastic semiclassical gravity ||, 0, ^, ||, similar to 
the expectation value of the stress-energy tensor in semiclassical gravity |Q. The noise kernel being a measure 
of the fluctuations of the stress tensor of quantum fields enters in a great variety of issues and problems ranging 
^ from the validity of semiclassical gravity to spacetime foams, from structure formation in the early universe to 
i^i fluctuations of the black hole horizon and trans-Planckian physics (for a review of this subject see ). Noise kernel 
^ ■ in hot flat space, one of the two examples considered here, has been studied by Campos and Hu Q]. It is useful 
7r| for performing a nonequilibrium quantum field theoretical analysis of the black hole nucleation problems beyond 
J> : Euclidean thermodynamics |lj, |l^ . 

In Paper I p^ , we have derived a general expression for the noise kernel of a quantum scalar field in an arbitrary 
^ curved spacetime as products of covariant derivatives of the quantum field's Green function. It is finite when the noise 
' kernel is evaluated for distinct pairs of points (and non-null points for a massless field). We also showed explicitly 
the trace of the noise kernel vanishes, confirming there is no noise associated with the trace anomaly. This holds 
regardless of issues of regularization of the noise kernel. 

The noise kernel as a two point function of the stress energy tensor diverges as the pair of points are brought together, 
sharing the ubiquitous ultraviolet divergence present in ordinary (point-defined) quantum field theory. To calculate 
a regularized noise kernel, similar to what was done before for the simpler stress energy tensor, it is desirable to have 
at the start an analytic and finite regularized expression for the Green function. When one can carry out a mode 
decomposition of the invariant operator and find an analytic solution to the mode functions there are established ways 
to proceed. For such cases the quantum stress tensor and its fluctuations can be determined using some regularization 
scheme, from the simple normal ordering |jl^, ^ or smeared field ||l^ in Minkowski and Casimir spaces to the elegant 
^-function or the powerful point-separation methods, as demonstrated for the Einstein Universe or for the Casimir 
effect by the authors earlier |l^ and others jl^, ^ Unfortunately, there are many important geometries 
for which an exact analytic expression of the mode functions is not available, such as the Schwarzschild black hole 
spacetime. 
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In this series of papers we work with the covariant point separation regularization method |22, |2J] in its modified 
form ^ to derive a finite expression for the coincident hmit of the noise kernel. The expression derived in Paper 
I for the noise kernel of a scalar field is completely general and can be used with or without consideration of the 
renormalization of the Green function. Also, the result there for the coincident limit holds for all choices of the Green 
function and the metric provided that the Green function has a meaningful coincident limit. In this paper and the next 
one, we apply this formal procedure to specific spacetimes of physical interest. We do this by working with an analytic 
form of the Green function. When such a form is available one can carry out an end point expansion displaying the 
ultraviolet divergence. Subtraction of the Hadamard ansatz expressed as a series expansion will render this Green 
function finite in the coincident limit. With this, one can calculate the noise kernel for a variety of spacetimes. 

An analytic form is obtained by invoking the Gaussian approximation introduced by Bekenstein and Parker [ p7| . 
For a massless scalar field in an ultra-static spacetime whose metric has an optical form (one where the Euclidean time 
T-time component of the metric grr = 1 ) this provides a closed expression for the Green function. In this paper we use 
the Gaussian approximation for the Green function for such quantum fields to evaluate the noise kernel in the following 
optical metrics : hot flat space, and the optical Schwarzschild spacetime, which is conformal to the Schwarzschild 
metric. For hot flat space, the Gaussian Green function is exact. For optical Schwarzschild, the Gaussian Green 
function is known to be a fairly good approximation for calculating the stress tensor p8| , which involves second 
covariant derivatives of the Green function. We will carry out this calculation for the noise kernel which requires up 
to four covariant derivatives of the Green function. Thus the validity of the Gaussian approximation will be tested to 
its new limit. A reliable check is provided by the trace of the noise kernel, which for massless conformal fields should 
be zero. 

Thus the goal of this paper is threefold: First, to present the detailed steps in the calculation of the regularized noise 
kernel for a quantum scalar field in a general curved spacetime using the modified point-separation scheme. Second, 
to derive the regularized noise kernel for a thermal field in flat space. Third, to determine the range of validity of the 
Gaussian (Bekenstein-Parker-Page) approximation to the Green function by examining the error in the noise kernel 
expression for the optical Schwarzschild spacetime. 

We present an outline of the calculation as follows: In Section f| we give a brief description of the Gaussian 
approximation to the Green function pTf for ultrastatic spacetimes . In Section [II we consider the regularization 
of the heat kernel in the class of optical metrics. We expand this Green function in an end point series so that it can be 
separated into a divergent piece and a finite remainder. The divergent term is independent of the approximation since 
we know this structure must be of a general form given by the Hadamard ansatz. After this is subtracted, we have 
a series expansion of the renormalized Green function. We then substitute this expansion in the general expression 
obtained in Paper I for the coincident limit of the noise kernel. The resulting expression is quite long and formal. At 
this point one can introduce the specific metric of interest and determine the component values of the Green function 
expansion tensors (by symbolic computation). From this we can readily generate all the needed component values 
of the coincident limits of the covariant derivatives of the Green function, along with the covariant derivatives of the 
coincident limits. These explicitly evaluated tensors are then substituted in the general expression obtained in Paper 
I to get the final result. We give two examples in Section For the case of hot flat space we derive the variance 
of the energy and pressure density for a quantum field at finite temperature. This is a useful compendium to the 
results obtained earlier for quantum fields in Minkowski and Casimir geometries in reference to issues like the 
validity of semiclassical gravity jl^, |l^, |2^, For a massless, conformally coupled field in the optical Schwarzschild 
(the ultrastatic spacetime conformal to the Schwarzschild black hole), we obtain for the regularized noise kernel at 
spatial infinity the same result as a thermal field in flat space, as it should, and a finite result at the horizon in a 
state conformally-related to the Hartle-Hawking state. However the latter expression computed with the Gaussian 
approximation has a non- vanishing trace. In Section ^ we study the nature and source of this error by examining the 
validity of the Gaussian approximation at successive orders. It works reasonably well to the third covariant derivative 
order. The inadequacy of the Gaussian approximation to the Green function for the calculation of the noise kernel 
arises from the Green function's failure to satisfy the field equation of the scalar field at the fourth order. 

II. GAUSSIAN APPROXIMATION 

We give a brief description of the Gaussian approximation to the Green function for quantum fields in optical 



spacetimes a la Bekenstein-Parker|27 and Page 



In the Schwinger-DeWitt proper-time formalism pi] the Green function is expressed in terms of the heat kernel 
K{x, 2/, s) via 



/>oo 

G{x,y) = / K 
Jo 



{x,y,s)ds. (2.1) 
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where the heat kernel satisfies 



d_ 

ds 



□ - 



The optical metric for an ultrastatic spacetime has the product form 



ds^ = gabdx°-dx = dr'^ 



9tj 



dx'-dx^ 



(2.2) 



(2.3) 



We assume in the Euclidean sector the imaginary time dimension is periodic with period 2tt/k — 1/T with T the 
temperature. For a black hole, n is the surface gravity but can be regarded as a temperature parameter here. This 
form of the metric allows the kernel to take on the product form 



with each of the kernels satisfying 



K(x,y,s) = Ki{t, T',s)K3{x,y,s) 



ir-l^A Ki{t,t',s) 



ds dr^ 



A'3(x,y,») = 



Equation ( p. 5a ) has the periodic solution 



Ki{t,t',s) ^ ^ ^ exp(-K 



^n^s + innAr) 



(2.4) 

(2.5a) 
(2.5b) 

(2.6) 



{At = t — t'). Equation (2.5b) in general is difficult to solve, but Bekenstein and Parker p7| find an approximate 
solution using the Gaussian approximation to the path integral representation. For K^, it takes the form 



(3)^1/2 

^i^3Gauss(x, y, s) = 3- exp 

(47rs)2 



(3^ 



(2.7) 



where ^^■'cr is the world function for the three-dimensional spatial geometry. In the above (•^W^^ is the VanVleck- 
Morette determinant for the spatial geometry. Since we have an optical metric, the four-dimensional world function 
is (T = '^^V -I- Ar^/2 and there is no difference between the three and four dimensional A-^^. 
Since the complete Hadamard-Minakshisundaram-Pleijel-DeWitt expansion [E2| would be 



n=0 

the Gaussian appr oxima tion is equivale nt to only taking the first term in this pow er series. 
By putting (2^) and (2/7) back into (2^) and carrying out the integration {'2A), Page obtains 



(2., 



GGauss(a;,y) 



sinh / 



Stt^t (cosh Kr — cos kt) 
as the Gaussian approximation for the Green function, where 



(2.9) 



III. NOISE KERNEL IN OPTICAL SPACETIMES 



In this section, the noise kernel for an ultrastatic spacetime with an optical metric is determined. For this class of 
geometries, we start directly with ( [2 .91 ) for the Green function. The first step is to expand this Green function about 
the coincident limit. Since the noise kernel has terms with at most four covariant derivatives, this expansion needs to 
be to fourth order in cr°, or second order in cr = aPap/2, and fourth order in At. Doing this expansion yields 

^ AW k2 , . 2 x 



4 



3^ (4 Ar^ - 6 Ar V + a^) I + O (ai , Jr^) 



By subtracting from the Gaussian Green function the Hadamard ansatz 



S{x,y) 



2Ai 



16 7r2 



(3.1) 



(3.2) 



(the V{x, x') term is absent since there is no log cr divergence present in the expansion of the Gaussian approximation 
to the Green function) we get the renormahzed Green function 



S. 



(3.3) 



The divergent term present (3.1) is canceUed by the divergent term from the Hadamard ansatz. 

Since our main interest here is to determine the coincident hmit of the noise kernel, we next turn to developing the 
series expansion 



Gr 



{An) 



-a^P a^^ a'^ G^p.rs 



(3.4) 



of the regularized Green function. With this, it will be straightforward to compute the coincident limits of the various 
covariant derivatives needed. We start by assuming the expansions 



A^ w l + a'Pa'«A(2)pg + a;Pa'«a'''A(3)„„^ + cr'P(j'«fT'''(T''^ A^") 



pqr 



pqrs 



At' 

Wi 

^ (0) 
W2 ~ W2 



a-P (T'« Sr^^^q + a'P a-"- Sr^^^gr + a'P a-^ a'"' a-' Sr^^^grs 



(3.5a) 
(3.5b) 
(3.5c) 
(3.5d) 



The specific values of the expansion tensors in these series are derived in the Ap pend ices. Carrying out the subtraction 
(|3.3D and substituting the expansions (3.5), we find the expansions tensors in (3.4) are 



g(0) ^ ^ 
^rcn 2 

G^i) - 

Gf(3) , - ^a(3) , + — At(3) . - io , u;^^^ 

'-'renofcc — g ^ abc ^ 45 2"° 1 

GitUcd = !^AW„,,, + ^(4A(2)„,5r(2),, + 4<5r(4),,,,-A(2),,5a. 
+ ^ (l6(Sr(2)„,^r(2),^ - 12 Sr^^^Ka gab + 9abgcd) 
-^Qab (w!^^ gcd + 2 w[^Kdj 
Using the explicit forms of the expansion tensor values, (|C8|), ( plCl| ), ( pl2a ) and ( |D5D , we get 



(3.6a) 
(3.6b) 
(3.6c) 

(3.6d) 



^rcndo 



72 

,2 



Rab 



45 



(3.6e) 

(3.7a) 
(3.7b) 



Gronabcd — 4320 Rab-cd + 5 Rab Red + 4 Rpaqb R/ d^) 

(12 T\b T\d - 16 T\b,c Sd^ + 4 Sa-" Sb^ Red - gab Red) 



2160 
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+37^ (16 Sa"^ Sb" 5 J 5/ - 12 5 J db" gcd + gab gcd) 



7560 
1 

'4 



^gab {w^2^ gcd + 2 wf^cd) 



where we have also used 



along with wj^^'a = 0, which hold for ultrastatic metrics. 
For flat space, these tensors reduce to 



(0) 



-240 (^-^ 



Rpq 



Rpqrs ) — 0, 



(3.7c) 



(3.8) 



^rcnao 



K 

180 


,6 



{-Vab + 'i'na'^ Vb'^) 



G^^labcd = 7^-::7:{'nab11cd-'i^'2ria'^Vb'^'ncd + iQr]a''r]b'^Vc'^Vd^) 

7550 



(3.9a) 
(3.9b) 
(3.9c) 



Now that we know the end point series expansion (3.4) of Grcn, t he c oincident limit of terms with up to four 
covariant derivatives arc computed. We simply differentiate the series (3.4) and then use the results from Appendix 
(^) for the coincident limits of the covariant derivatives of the world function a. The results are 



16 TT^ [G,en] 
16 TT^ [G„„;a] 
16 TT^ [G 

rcn;aoJ 
16 TT^ [Grcn;abc] 

16 TT^ [G^cn;abcd] 



q(0) 
rcn 

r,(0) 

^ rcn; a 

g(°) . « + 2 g(2) / 

^ron;(ao) ^ ^ron(ao) 

f- ( ^(2) I ^(3) 

" I ""rcn(afc:c) i ^rcn(ofcc) 



6c 



(3.10a) 
(3.10b) 
(3.10c) 

(3.10d) 



12 



(2 G^cn(Qbc;d) + G^^^(afc;c<i) + 2 G^^^(ofccd) 



G^-^npa {Rbcd^ — 2Rbdc') + G^^l^pb {Racd^ — 2 Rad/) 
+ G'rcnpc {Rabd^ ~ 2 Radb^) + Gj-.^^jprf {Rabc^ — 2 Racb^)^ 
~^G^f^^-^abcd 



(3.10e) 



We now have all the information we need to compute the coincident limit of the noise kernel (see Eqn (3.24) of 
Paper I). Since the point separated noise kernel Nabc'd'{x,y) involves covariant derivatives at the two points at which 
it has support, when we take the coincident limit we can use Synge's theorem to move the derivatives acting at the 
second point y to ones acting at the first point x. Due to the length of the expression for the noise kernel, we will 
here give an example of the calculation by examining a single term. The complete expression for the coincident limit 
of the point separated noise kernel can be found in Paper I as Eqn (4.16). 

Consider a typical term from the noise kernel functional (Eqn (3.24) of Paper I): 



Grcn ; d b Gi cn \d' a ~^ Gi cn ; c' a G^cn -^d' b 

As was derived in Paper I, the noise kernel itself is related to the noise kernel functional via 

Nabc'd- = Nabc'd' [G,.e„(x, y)] + Nabc'd' [Greniy, x)] . 



(3.11) 



(3.12) 



We account for this by adding the same term but with the roles of x and y reversed. Taken together, we need to 
analyze 



Grcn; c' b G^cn; d' a ^ Gi-cn;a' d G^-cn; 6' c ^" G^cn ; c' a Gj-cn; d' b ^" G^cn; a' c Grcn \b' di 



in particular, its coincident limit: 

[Grcn;c'6] [Grcn;d'a] [Grcn;a'd] [Grcn;6'c] [Grcn;c'a ] [Grcn;ti'fc] + [G 

We apply Synge's theorem to remove any explicit reference to the point y, 

([Grcn:a]:d ~ [Gi-cn:ad]) ([Grcn:6];c ~ [Grcn:6c]) 



[G 



ren:b' d\ 



(3.13) 



(3.14) 
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+ ([Grcii:d];a ~ [Gron;ad]) ([Grcn:c]:6 ~ [Grcn;6c]) 
+ {[G ren;aj ; c L^rcn ; a 

+ {[Gien;c];a — [Grcn-ac]) {[Gien;d\:b — [Grcn;bd]) (3.15) 

and express the results in terms of the expansion tensors ( |3.10e| ) and find 

8 (^G[:J^ad G^J^bc + G[:J^ac G^J^bdj (3.16) 



Though this may look relatively compact, when we use the results (3.7c) to get the final form, in terms of the local 
geometry, this is no longer the case. Making these substitutions, our pair of terms become 

{Rad Rbc + Rac Rbd) 



2592 7r2 

^ { (4 (Jfc"" (5d^ - Sbd) i?ac + (4 (5b^ 5 J ~gbc) Rad 



12960 7r2 

+ (4 6a^ 6/ - gad) Rbc + (4 <5a^ SJ - gac) Rbd} 
+ a.onn 7 Sa'^ Sb' 5/ 5d^ + gad gbc + gac gbd 

64800 TT^ 

-4 {5b"' 5/ gac + h'" 5 J gad + 5a'' 5/ gbc + <5a^ 5 J gbd)} (3.17) 

The noise kernel functional consists of 25 such terms, some quite a bit more involved. This is especially true when a 
single Green function has four derivatives acting on it. To gain insight into the physics we work with some specific 
spacetime. Q 

When we do choose a metric, we proceed by directly evaluating the components of the expansion tensors (3.7c). 
Once these expansion tensor components are known, the actual components of the coincident limits of the covariant 
derivatives of Grcn are in turn computed using (B.lOe). Then it is straightforward from there to get the covariant 
derivatives of the coincident limits, since application of Synge's theorem will move derivatives acting at y such that 
we also need the covariant derivatives of the coincident limits. Now that we have the component values of all the 
covariant derivatives of the various coincident limits of the regularized Green function Grcn, we substitute them in 
the coincident limit expression for the noise kernel, Eqn (4.16) of Paper I, and arrive at the final result we seek. 

Before turning to specific metrics and as a check, we can reproduce the derivation of the rcnormalized stress tensor. 
We start with the point separated expression for the stress tensor, which for a massless, conformally coupled scalar 
field is 

{Tab{x, y)) = i (^gP b Grcn-p'a + g'' a Groiup'b) - ^.9^ q Grcn;p''' 9ab 



6 
1 

^6"'"' V 2 

We take the coincident limit and utilize Synge's theorem to obtain 



^G^enf ^^J^^'K Rab] (3.18) 



(T'afc)i.on ~ g [Grcn:a];b + 3 [Grcn:6];a ~ [Groii];a6 ~ 6 [Grcn;afc]) 
7i ('^ [Grcn;??] '^ [Gi-cn];p^ 3 [Grcnip^]) gab 
— Y^[Gren] {Rgab — Rab) 



(g(°) R - 2 gW - 12 G<,ll/) Sab (3.19) 



With the explicit values ( 3.7c ) for the expansion tensors we recover the familiar result: 



^4 



(Tab),,,, = {gab - 4 Sb^ (3.20) 
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IV. EXAMPLES 
A. Hot Flat Space 



The first example we consider is that of a finite temperature T = k/2tt quantum scalar field in fiat space. With 
this, the stress tensor takes the usual form 



/ , , ,111, TT^T* 

(r.,)=diag{--,--,--,l}p, P = -^ 



with a;" = (x, y, z, t) 



Using (3.9c) and (3.10e), the non-zero coincident limits of the derivatives of Gren are 

T2 



[Gron] — 



12 



[G^cn-ab] = -TZ {5ab~^5J6b^) 



[G 



rcn ;afccdj 



90 
315 



{Sab Scd - 12 (5a^ Sb'^ Scd + 16 Sa'^ 5b'' 5 J 5/) 



(4.1) 

(4.2a) 
(4.2b) 
(4.2c) 



Taking the coincident limit of the massless case for the noise kernel, Eqn (4.16) of Paper I, (for this case, we keep ^ 
arbitrary), we find 

Nabcd = ^^{32 (7- 28^+162^2) g-a 3^.9%. 

-8 (7-42^ + 130^2) {gcdg\9^b + gabg\ fd) 

-4 (7 - 28 C + 148 ^2) (gf,^ g\ g^^ + g^ ^ g^ ^ + g^^ g\ g^ ^ + g^^ g\ g^ ^) 

+4 (7 - 63 e + 167 f ) gab gcd + (7 - 28 ^ + 108 C^) (g^^ g^^ + g^^ g^^) } (4.3) 

From this, we can immediately compute the trace: 



(4.4) 



1350 

which we see vanishes for the conformal coupling case (^ = 1/6)- The non- vanishing components for general coupling 
are 

_4 rT,S 



^ TTTT 

^ xxxx 
^ xxyy 

-^TXTX 

^xyxy 



1^(7-14^ + 270^^; 
'21 - 154^ + 442^2) 



113400 

^4 rpS 



56700 



56700 
226800 

^4 j,8 



(7-2U + 93C2) 
(7-63^ + 167^2) 

(21-84^ + 484^2 
(7-28^+108^2) 





^4 j,S 




5400 
7r J 


5->o 


5400 




8100 

TT J 




8100 
7r _/ 




10800 

^4 rp8 


32400 



R {Trr) 



± (T ) 

2 XX I 



— 3 {^tt} {T'xx) 
= l(Txxf 

4 {Ttt) {Txx) 

- - (T \2 



(4.5) 



226800 

and those that follow from the symmetry of the metric. 

We use the dimensionless measure of fluctuations |l^, ^ |l^; (this is not a tensor, but a measure of the fluctuations 
for each component): 



^abcd 



{TabTcd) — {Tab) {Ted) 



bed 



4^Nabcd + {Tab) {Ted) 



(4.6) 



{TabTcd) 

From inspection, < Aabed < 1- Only for A ~ can the fluctuations be viewed as small. For A ~ 1 the fluctuations 
are comparable to the mean value. 

For hot flat space, the results for A are 



abed Aabed 

2 (7-145+270^^) 
35-28 4+540 4^ 
2 (21-1544+4424^) 
XXXX 49-308 4+884 4^ 
4 (7-214+934^) 
TTXX 49-84 4+372 4^ 
4 (7-63 4+167 4^) 

^■'^yy 35-2524+6684^ 





= 




1 

~ 6 


2 
5 


- 0.4 


73 

136 


- 0.54 


6 
7 


^ 0.86 


137 
200 


- 0.69 


4 
7 


^ 0.57 


73 
136 


- 0.54 


4 ^ 
5 


- 0.8 


41 
104 


- 0.39 



(4.7) 
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From these results we see, even for the simple case of thermal fluctuations in flat space, the fluctuations present in 
the stress tensor are important. Discussions on the implication of the fluctuation to mean ratio can be found in 




B. Optical Schwarzschild black hole 

We now consider the optical spacetime conformally related to the Schwarzschild black hole spacetime. For this 
spacetime, the line element is 

ds^^dT^+(^l~^^^ \r^+(^l-^i^ \^d0^+sinHd(l)^) (4.8) 

Taking k = 2ttT and T — 1/{8ttM) we choose the quantum state corresponding to the Hartle-Hawking state in the 
conformally-related Schwarzschild spacetime. We use the spacetime coordinates x° = (r, 9, 0, r) and introduce the 
rescaled inverse radial coordinate x = 2M/r = l/(47rTr). Spatial inflnity corresponds to x = and x = 1 is the black 
hole horizon. For a massless conformal scalar fleld {m — 0, £^ — 1/6) the stress tensor is 

(T/)=diag|-,-,-,-l|^ (4.9) 

We recover the standard thermal result for the stress tensor. The component values for the noise kernel are 

Nr^r^ = (657 - 1050 xV 8400 x^- 16800 x^- 16997400 a;** 

2041200 ^ 

+80206000 x^ - 140910000 x^" + 109242000 x" - 31515750 x^^) (4.10a) 

Nr\'' = — ——r (137 + 1120x3 -210x^^ + 80640 x^- 103600 x^ + 16800 x^ 
2041200 ^ 

+9829800 x^ - 43285200 x^ + 73602000x^0 - 57186000 x" 

+17057250x^2) (4.10b) 

iV/e^ = ^^^^^^ (137 - 560 x^ + 1470 x'' + 30240x5- 44800 x^- 16800 x^ 
-15805800 X* + 73203600 x^ - 127386000 x^" + 98820000 x" 

-28788750x^2) (4.10c) 

iV^^/ = (-657 + 6720 x=*- 5670 x''- 42000 x^ + 484920 x® + 4202480 x^ 

6123600 ^ 

-12514800x^0 + 9439200 x" - 1548450x^2^ (4.10d) 

7^4 rp8 

^r^/ = ^ ^o^^^ (-657 - 3360 x^ + 4410 x^ + 8400 x^ + 25200 x^ - 53478360 x^ 
6123600 ^ 

+248828960 x^ - 434589600 x^° + 337051800 x" - 97825050 x^^) (4.10e) 

_-4 rp8 

N/e^ = (123 - 5040 x^ + 3150 x* - 120960 x^ + 176400 x^ - 25200 x^ 

6123600 ^ 

+2858760 X® - 10502560 x^ + 15885600 x^° - 12447000 x" 

+4178250x^2) (4.10f) 
From the component values of the noise kernel we can compute its trace 

4 tt"* x^ 

N = N/a'' = (9720 - 45832 x + 80520 x2 - 62424 x^ + 18009 x^) (4.11) 

567 

We know from prior results that this should vanish, since we have worked with a massless conformally coupled scalar 
fleld. Thi s fa ilure of the trace of the noise kernel to vanish is due to the failure of the Gaussian approximated Green 
function (2.9) to satisfy the fleld equation to fourth order. 



V. FAILURE OF GAUSSIAN APPROXIMATION AT THE FOURTH ORDER 



To be sure that this error does not arise from the symbolic manipulation, let us mention ways to check the correctness 
of the algorithm. The basic procedures for generating the needed series expansions are recursive on the expansion 
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order. (The recursion formulas of the expansion used in point-separation are collected in the Appendices). For the 
noise kernel we need results up to fourth order in the separation distance. The well established work for the stress 
tensor is to second order. This provides a check of our code by verifying we always get the known results for the stress 
tensor expectation value. Once we know the second order recursion is correct, we know the algorithm is functioning 
as desired. The correctness in the (new) fourth order terms becomes particularly important in the work in Paper III, 
when we consider metrics conformally related, as we get intermediate results of up to 1100 terms in length. 

We can check the accuracy of the Gaussian approximation by using the computed component values of coincident 
limit of the covariant derivatives of Gren- We have assumed 



Gr, 



Gicn R 

6 







(5.1) 



This can be used to test the approximation order by order. To test up to the second order, we just take the coincident 
limit 



R [Grc 



= 0. 



(5.2) 



For the metric (4^), the scalar curvature is 



R = 



6M^ 



= -24 TT^ x'^ 



(5.3) 



while the results of our computation of the component values of the coincident limit of the covariant derivatives of 
Grcn yield 



4 TT^ 16 TT* x"^ 

[Grcn] = 5 , and [Grcn;p^] = 5 ■ 



(5.4) 



With these values, Eq. ( f3.2| ) can be seen to be satisfied. Thus the Gaussian approximation is good to the second order. 
(This had better be the case, as the approximation at this order has been checked against numerical computations of 
the stress tensor by other authors. See e.g., description in js^, |33| ) 

To check the third order term, we take one covariant derivative and then the coincident limit of (fe.ll): 



Using the results [Grcn;a] = 0, 



[GrGn;p''a] ^ g {R-.a [Grcn] + R [Grcn;a]) — 



(5.5) 



R.^ = {384 7r^T3a;^ 0,0,0} 



(5.6) 



and 



p /• 



^5 7^5 ^5 



-,0,0,0 



(5.7) 



we see Eq. ( p.5[ ) is also satisfied. This has to be the case: (3.7c) shows the third order expansion tensor does not 
have any contribution from the W{x^y) part of the Hadamard ansa tz a nd this is the only place a lack of symmetry 
in Grcn (a;,?/) could come in. Therefore the Gaussian Green function ( ^.9[ ) is symmetric. For symmetric functions, the 
odd order expansion tensors ar e de termined completely by the even order tensors (see (B16a)). 
Continuing to fourth order, (5.1) becomes 



[Gre„;p^a''] " 7 {R;a' [Grcn] + -R''' [G 

rcn ; a ] + R,a [G ICll 

+ R [Grcn:a'']) = 
[Grcn;p''a''] ^ g {,R,a' [Grcn] + R [Grcn;a'']) = 



(5., 



Proceeding as before and evaluating the left hand side above, the component values are 
diag{ - 128 .«T«.'^ 



315 (162648 

128 7r° T° X 



315 



746888 X - 
(8424- 29704 a; 



1295880 x^ - 1005480 x^ + 293805 x^) 
44040 x^ - 34560 x^ + 11835 x'^) , 
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(8424 - 29704 x + 44040 - 34560 + 11835 a;'* 

128 tt" T° x" (aynn _ /| Koon „ i onc^On -,-2 _ ROAOA ^'i J_ 1 SinnQ -,,4 



315 



189 



9720 - 45832 x + 80520 - 62424 a;^ + 18009 x^) } (5.9) 



The failure of the left hand side of ( |5.8D to vanish shows that the failure of the trace to vanish comes from the 
limitations of the Gaussian approximation. The Gaussian approximation is only useful up to the third order in cr". 



With this knowledge, the trace (4.11) becomes our measure of the error in the noise kernel from the use of the 



Gaussian approximation. It is important to note that the noise kernel trace N vanishes as a; — > 0, or, r — > oo, i.e. 



where one would expect the effects of curvature to vanish. We can also see from ( 4.101 ) that the noise is finite at the 
horizon [x — 1). 

Using our derived expression for the no ise kernel we see that its trace vanishes at spatial infinity, thus we can trust 
our results there. Using the measure (^^), the fluctuations at r ^ cxi are 



abed : tttt rrrr 6969 rrrr tt69 rr6 

73 137 137 J73_ 73_ 41 
136 200 200 136 136 104 



A . J73^ 137 137 73_ 73_ 41^ (5.10) 

^abcd 



which match exactly the results (4.7) for hot flat space with conformal coupling, another reassuring fact. Since the 



computation of the noise kernel for metric (4.8 ) is much more involved, this provides yet another check of our symbolic 
computer code. 

Now having shown that it is truly the Gaussian approximation that is a t faul t for the failure of the noise kernel 



trace to vanish, we can use its value as a measure of the error in the results ( 4.10i ). Since N should be zero, N/Nabcd 
is a dimensionless measure of this error. At the horizon (x = 1), 

abed : tttt rrrr 6969 TTrr tt96 rr69 

^44 : 627% 791% 791% 1390% 1390% 19855% ^^'^^^ 

These errors show the Gaussian approximation fails to provide reliable results for the noise kernel near the horizon of 
the optical Schwarzschild metric. We expect this be the case also for the Schwarzschild metric near the horizon, an 
explicit calculation will appear in the next paper. In the above we have identified the occurrence of significant error 
begins at the fourth covariant derivative order. 

In conclusion, towards the three goals set for this paper, we have detailed the steps in implementing the modified 
point separation scheme under the Gaussian approximation for the Green function for a quantum scalar field in the 
optical spacetimes. We have derived the regularized noise kernel for a thermal field in flat space. We have obtained a 
finite expression for the noise kernel at the horizon of an optical Schwarzschild spacetime and recovered the hot flat 
space result at infinity. From the error in the noise kernel at the horizon we showed that the Gaussian approximation 
scheme of Bekenstein-Parker-Page applied to the Green function which provides surprisingly good results for the stress 
tensor involving the second covariant derivative order of the Green function, fails at the fourth covariant derivative 
order. 
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APPENDIX A: WORLD FUNCTION a 



In this appendix, we review the properties of the world function a. Wc also demonstrate how symbolic computations 
are implemented and used in this work. Christensen's method for determining the coincident limit of covariant 
derivatives of functions defined via a covariant differential equation is reviewed. 

The world function is defined to be one half of the square of the geodesic distance between two points on a differential 
manifold. As such, it satisfies the equation 

-^=0 (Al) 

along with the initial value 

[fT;a]=0 (A2) 



To determine [ciab], we take two derivatives of (Al). To help illustrate how these calculations are done on the 
computer, the output presented here is direct output from MathTensor [p4| , with Mathematical^ carrying out the 
formatting. We use the MathTensor function CD and get 

CT;ab ^ ^ ^ ^ - U {A6) 

This is put into a canonical form via Canonicalize 

(T-ab- cr-pacr-h'' - o-;pa.J'i, = ^ (A4) 



The condition (A2) is encoded into Mathematica by defining a rule 

Ci [ CD [sigma, laj ] -> 0, where Ci[] is a function defined to represent the coincident limit and is formatted to be 
displayed using the standard [• • • ] notation. Using this rule, the coincident limit is 

[a;afc] - [cT,pa] [<J,b'] = (A5) 

This immediately shows [a;ab] — gab and such a rule is defined. 
Proceeding with one more covariant derivative, 

~ Cr-abc — <^-pc<^-J'b — 0';ph(T;Q^c — <7-pa<7;b''c " CT -p <7;a^ be, (A6) 

and using the two rules we already have, we recursively get 

= [a.abc] + [cr.acb] (A7) 

Using MathTensors' OrderCD, which commutes the covariant derivatives on each term until they are in alphabetical 
order, we get the result 

- 2[a,abc] + [cr;p]Ra^bc 

= 2[a,abc] (A8) 

where ( |A2|) is used to go from the first to the second line. We have the coincident limit of three covariant derivatives 
acting on a vanishes. 

Proceeding with one more covariant derivative: 

= 0':padO':b^c + Cr-pacO':b^d + Cr-pabCr-c^d — Cr:abcd + '^■,pd'^;a^bc 
+0" :pc O- -a^bd + -pb Cr -a^cd + ipa O' -b^cd + <^;p (^-a^bcd 

= [cr.abcd] + [j-acbd] + [(7;adbc] (A9) 

where once again we use the rules we already know. Now commuting the covariant derivatives 

= 'i[(y-abcd] + {Ra'bc-d + Ra'bd-c) [(^-.p] 

+ [CT;pd] RJ'bc + [fT;pc] RJbd + [<J-pb] Ra^cd + [<J;pa] Rb^ cd (AlO) 

and using the known rules, we get the equation 

= 3 [(T;a6cd] + Racbd + Radbc (All) 
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which we can solve for the term we need. In practice, Mathematica's function Solve [] is used, giving us the rule 

[f ;abc<i] ^ {Racbd + Radbc) (A-12) 

By knowing the coincident limit of n — 1 covariant derivatives of cr, we determine the coincident limit of n covariant 
derivatives. This is the recursive algorithm developed by Christensen. It is the main idea we use for computing the 
expansions needed in this work. A general outline starts by assuming we have the rules for n — 1 covariant derivatives, 
then, 



1. Take n covariant derivatives of the defining equation (in this case, Eq (Al)); 

2. Use the rules for n — 1 covariant derivatives to get the coincident limit; 

3. Commute the covariant derivatives; 

4. Use again the rules for n — 1 covariant derivatives on the terms generated; 

5. Solve for the coincident limit of the n derivative term; 

6. Define a new rule for this term. 

These steps are iterated until all terms needed are generated. 

For the world function a, we need to carry this out to eight covariant derivatives. The seventh and eighth order 
(derivative) terms become quite large. In fact, when computing these expression, we only substitute (steps 2 and 4 
above) up to four covariant derivatives and still get results with 240 terms for the seven derivative result and 1101 
for the eight. We only finish carrying out the recursion when we use these highest order values. 

The results for five and six covariant derivatives are: 

[c^-abcde] = — ^ {Racbd-e + Racbe;d + Radbc-e + Radbe;c + Raebc-d + Raebd;c) (^13) 

[c^iabcdef] — ~~ {Racbd:ef + Racbe;df + Racbf [de + Radbc;ef + Radbe;cf 
+ Radbf;ce + Raebc-df + Raebd-cf + Raebf-.cd + Rafbc-de 
+ Rafbd-ce + Rafbe-cd) /5 

{.Rpecd Rabf^ ~^ Rpebd Racf^ ~^ Rpebc Radf^ ~^ Rpead Rbcf^ 



y^^pecd ^^abj ' ^*^pebd ^^acf ' -'-^^pebc -'-^adf ' -'-^^pead -'-^^bcf 
+ Rpeac Rbdf^ + Rpeab Rcdf^) /9 
+ {Rpdcf Rabe^ + Rpdce Rabf^ + Rpdbf Race^ + Rpdbe Racf 
+ Rpdaf Rbce^ + Rpdae Rbcf^) / 45 

+ {Rpcbf {Rade' + T Raed^) + Rpcbe [Radf^ + 7 Rafd"^ 
+Rpcbd {Raef^ + 7 Rafe^) + Rpcaf {Rbde^ + 7 Rbed^) 



(A14) 



APPENDIX B: END POINT SERIES EXPANSION 



The basic input into the computation of the stress tensor or noise kernel is the Green function, a perfect example 
of a bi-scalar. We want to express it in such a way we can easily identify how it depends on the distance between its 
two support points. This leads us to consider series expansions of bi-scalars. The techniques are also useful for the 
series expansions of bi-tensors. 

The world function a introduced above provides the ideal geometric object for such a construction. It contains 
both distance and direction information. For a bi-scalar S{x,y), the end point expansion is 

Six, y) = + (jP^i) + + . . . + • • • (jP-^^-p. + ' ' ' > (Bl) 

so-called because the expansion tensors A^a^...a^ = ylQ"^..^^^ (a;) have support at one of the end points for which S(x, y) 
has support. 
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It is only the symmetric part of the expansion tensors A^\.a„ that contributes to the expansion, since they are 
contracted against symmetric products of cr^* 's. Moreover, the expansion tensors are order n in distance contribution 
to the bi-scalar S{x,y). We also find it convenient to have an expansion where the distance dependence is separated 

from the direction dependence. To this end, if is the unit vector along the geodesic from x to y, cr'^ = {2a)^ p", 
the expansion (Bl) can be re-expressed as 



S{x, y) = ^(0) + CTi^(i) + (tA(2) + • • ■ + (jt + • . 



(B2) 



where A^") = 2tpPi . . ■pP^A'p^'...p^. Now the expansion scalars carry the direction information. 

When multiplying series, this form readily collects terms by their order in distance. In the context of symbolic 
manipulation of series on the computer, this alternate form greatly improves processing speed. 

The first series expansion to consider is that for the world function, which, by virtue of its defining differential 
equation is given by 



a{x,y) 



9pq 



(B3) 



This is exact and from it we can see that all expansion tensors of order n > 3 vanish. This in turn tells us that 
[''';(aia2a3 -a„)] — 0, the Coincident limit of three or more symmetrized covariant derivatives of the world function 
vanish. This can also be seen by direct inspection of the previously given expressions for the coincident limits. 

We now turn to relating the expansion tensors to the coincident li mits of the derivatives of the scalar S. This is 
done by taking covariant derivatives and then the coincident limit of (Bl). We immediately get 



Ad 



Taking one derivative and the coincident limit gives 
Two covariant derivatives yields 

We only need the symmetric part of A^^^i,, or if wc assume A'^^^afc is symmetric, we solve for: 



[S]:ab + ['5':ab]) 



We write this as 



^ ab 



-[S-a]-b 



[S];ab [S-ab] 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



where we use the standard notation = to denote equality upon symmetrization. In terms of symbolic processing, this 
is implemented by taking each term of a tensorial expression and putting all free indices in lexigraphic order. We also 
define rules that set to zero any Riemann curvature tensor Rabcd when either the first or second pair of indices are 
free. For example, consider when we take three covariant derivatives: 

+A^^Kb;c + A(2)ac;. + A^^\,,, + A^^\r^,a + A^^K,,, + A^^K,,, 
+A('\-M + A^'\:ac + A^'K-^ab + ^^^Ijf^ + ^^^^ 

+A^°^;abc (B9) 

Now putting all free indices in lexigraphic order and then using the above rules for the Riemann curvature tensor 
gives 



[S;abc] = 6A^'\bc+6A^^\b;c + 3A^'\,bc + A''°^,abc + 
= 6 -A^ ^afoc + S-A'- ^a6;c4"3-A^ ^a:6c+A.^ ^ -.abc 



2 ^ p Rabc^ 



(BIO) 
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This can now be solved for A^^\i,c and the previously determined results for ^^"^ , A^'^^a and A^'^\b used to get A^^\bc 
solely in terms of the coincident limit of up to three covariant derivatives acting on 5*. Nothing new is encountered 
when determining the rest of the expansion tensors. We now give the results for up through A^^\bcdefgh- 



bc\ 



ab\ • c 



be 



6 



[S]:a 



be 



2^ [Seabed] 



g[>5':a6c];d + ^[>5';a6];cd ~ g[5';a];hc!i + ^['S'];abc(i 



abode 



1 , 1 



12 



bc\:de 



12 



[S;ab]- 



cde 



24 



[S-A 



bcde 



120 



abode 



bcdef 



'^[S-^a];bcdef 



48 
1 



abcdefg 



■,ab];cdef ^ 
;abcdef g] 



120' 
1 



4 ■.abcd]:e f ^^{S -abcl^de f 



(Blla) 



(Bllb) 



(Bllc) 



(Blld) 



1 



[S-^abcde]-J g 



'T7j['^'°'"='*]'^/9 ~^ T77['^;a''c];de/g ^ -^^[S ;ab];cde f g 



144 
1 

720 

1 



abcde f gh 



\a\ ]bcdefg 
-^abcde f gh\ 



144' 
1 

5040 
1 



[S] 

\abcde f g 
\S -^abcde f g\ -.h 



(Bile) 



1 



" \S -^abcde f \ ;gh 



40320' 5040' ' 1440' 

~':^^[S;abcde]-,f gh + ['5';abcd] ;e/g/i ~ l^^iS ■abc];def gh 

1 , 1 , 1 



1440' 



IS: 



ab\ \cdefgh 



5040 



bcde f gh 



40320 



abcde f gh 



(Bllf) 



These relations simplify considerable if the scalar S is symmetric, S{x,y) — S{y,x), for the symmetrized odd 
derivatives are determined by the even derivatives. For one derivative, 



S{x,y)-a' = S{y,x).a' ['S'la'] = [5';a], 

and applying Synge's theorem to the left-hand side above yields 

[S];a 



For three derivatives, we have 



Once again, using Synge's theorem and (B13): 

[<5';ab];c + [S-ac];b + [S-bc 



[S-a'b'c'] — [>5';a&c] 
[S]-acb 



[S-bca] — [>S';abc] 



It follows from this that 

4 [5';(afcc)] = 6 [S';(a6];c) — [S]-(abc) 

The results for five and seven derivatives are 

2 [S';(afcc<ie)] = 5[S-(abcd];e) — 5[S-(ab];cde} + [S]-(abcde) 

8[S.,^abcdefg)] = 28 [S ,^abcde f] .g) - 70 [S , (abed] .ef g) + SA [S :(abj \edef g) 
— 17 [S]-(fibcdefg) 



(B12) 

(B13) 
(B14) 

(B15) 

(B16a) 

(B16b) 

(B16c) 
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With these results, the equations for the even expansion tensors for a symmetric function simplify: 

= [S] 

4!AWafccd = [S.abcd] 
6lA^-^\bcdef = [S-abcdef] 
SlA'^^^bcdefgh = [S-abcdef gh] 

while the odd expansion tensors are given in terms of the even tensors: 

6\A^''\bc.de = -360 A^^\hcd;e+ 50 A^^Kb;cde- 3 A^°l^abcde 
S^-A'^^Kbcdefg = -20160 ^(6),,,,,,^.^ + 1680A^^Khcd;efg ~ 168 A^^\h-cde f g 
+ 17 A^^^ -abode f g 



Very often, we need the covariant derivative of a series expansion (Bl): 



pq-,a 



(B17a) 
(BlTb) 
(B17c) 
(B17d) 
(B17e) 

(B18a) 
(B18b) 
(B18c) 

(BlSd) 



(B19) 



If we replace a-ab with its series expansion, then we readily have the series expansion of S-^a- We can get this via the 
above relations by merely replacing S with a-^ab- In particular, we want the expansion 



+B'^^Kbpgrs a'" a'"- a'' + B^^'Kbpgrst a'P a'^ a''' a'' a'' 
I n(6) 

+£> abpqrstuCr Cr Cr CT' (T' (J 



We immediately have 



and 



For the second order term: 



B^°\b = [cr-,ab] = gab 
B'^^Kbc = -[(T-,ab]-,c + [a-;afcc] = 
^(2) _ r_ 1 I W;ab];cd , [o';Qbcd] 



ibcd 



g {Racbd + Radbc) 



(B20) 
(B21) 

(B22) 
(B23) 



Now we have to be careful about how we carry out the symmetrization: it is only the indices c, d that are contracted 
over in the series (B20). So it is only the free indices other than a and h in this and the following that we put in 
lexographic order (our routine for ordering free indices can be given a list of indices to exclude from ordering): 



r(2) 

i> abed — — 



Ra 



cbd 



(B24) 



Equality upon symmetrization of all indices but a, h is denoted by The rest of the expansion tensors are computed 
in the same way; the results are 

/ Racbd-^e 



^(3) 



5(4) 



abcde 



abcdef 



abcde f g 



5(6) 



12 

1 

180 
1 

360 
1 



abcde f g h 



15120 



( — 3 Raebd-.ef + 4 Rpead Rbef^) 

{Raebd;e f g 3 Rpdbc-.e Rag 3 Rpdac\e Rbgf^^ 

( — 219 Rpdac-,e Rbfg^;h — 6 Racbd-,efgh 

— 114 Rpdbc:ef Ragh^ — ii4:Rpdac:ef Rbgh^ 
+ 24 Rpeaf Rqhbg Rc^ d'' + 8 Rpeaf Rqhbg Rc"^ d^) 



(B25a) 
(B25b) 

(B25c) 

(B25d) 
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APPENDIX C: VANVLECK-MORETTE DETERMINANT 



Other than the world function, the other main geometric object we need is the VanVlette-Morette determinant, 
defined as 



D{x,y) = -det(-(T;afc' 
In the context of Green function, what appears is 



(CI) 



which we focus on. Using 2a = &''^(J-^p the VanVlette-Morette determinant is seen to satisfy 

D-^ {Da'P).p = 4 => (4 - Dcr) - 2AY/CT'f = 

along with 



[D] = 9{x) 



Ai/2 



from which we readily get 



0. 



(C2) 



(C3) 



(C4) 



(C5) 



We could at this point proceed as we did with a to determine the coincident limit expression of covariant derivatives 
of A^/^. But what wc need is the end point expansion of A^/^(a;,?/) to sixth order in cr". With this in mind, we set 
out to directly determine the series. We start by assuming the expansion 



Ai/2 = 1 + A(>V^ + A(3),cT^'aV + A(4),,aWV^- + A^^^^aVVVV* 



(C6) 



and substitute back into (|C^). We use (B20) and cr" = (Scr)^ p". The expansion tensors Ai"?..a„ are determined by 
collecting terms according to their order in a and setting to zero. 



From (|B20|) and (B25d), we have the series expansion 



u t\ 

s I 



^ -y PP P'l P"- P' (219i?p,,„;,i?/t'";„ - 6i?; 



1890 



-L9A F? F? ^ TD V X 

~r^^ J^pvqw sx -^t u 

J-S R F? ■"^ T) X v\ 

-rO J^pvqw sx -'^t u ) 



'pq;rstu 

99Q TD ID W V 



(C7) 



We have made the split into cr and p"" so we can readily carry out the needed multiplication in the series (C6), once 

1/2 I I 

we put it into the same form. The last piece we need is lS.-a ■ This is obtained by differentiating (C6) and then 



s ubs tituting (B20) for the a-ab terms that arise. Once this is done, we have all the terms for the series expansion of 
(C3). The order a term: 



2apPp'i 



(l2 A(2)p^ - 







A(2) 



Rah 
'12 



(C8) 



The order a 2 term starts off as 



(J 2 pP pi p^ 



3V2 



^24 ^ pq;r H~ Rpq;r H~ '^2 ^ pqr^ — 
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Using (|C8|) shows 



For the order term, we find 



a(3) ^ ^ 



24 



bed 



(-360 A(3)„fc,.<i + 3 i?afc;cd + 60 A^^),^ Red + 120 A^^)^, i?^. 



1440 



+ 120 A(2),p R,P^^ + 4 



g ,p 



The fourth and fifth terms above vanish since we only need equality up to symmetrization. Using (| 
the final form for this term becomes 

A^^'-'abcd = ^440 ^ab-cd + 5 Rab Red + 4: Rpaqb Rc"^ d^) 

The last two coefficients are determined in exactly the same manner. The results are: 

1 



a(5) , . 

^ aocae 
M'\bedef 



1440 
1 

362880 



(4 Rab-ede + 5 Rab-.c Rde + 4 Rpaqb-.c Rd''e^) 

(315 Rab:c Rde J + 180 Rab-.cdef + 378 Rab-ed Ref 
35 i?ah i?c(J Ref + 84 i?Qf, Rpeqd Re"^ — 270 Rpaqb;c Rd'^ e^ J 
^288 i?pa96;cd i?e«/P + 64 Rpaqb Rrde" Re'' f) 



(C9) 



(CIO) 



(Cll) 



and dClOl ), 
(C12a) 

(C12b) 
(C12c) 



APPENDIX D: SERIES EXPANSION FOR Ar 



In this section we determine the expansion 



Ar^ = CT'P (T'« (5T(2)p^ + CT'P CT''^ <5r(3)p^^ + a'" o"-'' o--" (Sr^^'^p,^,. 



(Dl) 



Since this is the expansion of a symmetric function, the expansion tensors are related the coincident limit of the 
covariant derivatives of At^ via (B17e) and (B18d). We also use [At] = 0, justifying the series expansion starts at 
order two in ct°. The expansion tensors are 



c (4) _ [Ar;^b] [AT;ed] , [Ar^g] [Ar.fc.rf] 

'^'^ abed — A ~^ 



(D2a) 
(D2b) 

(D2c) 



To evaluate the covariant derivatives, we start with Ar^a — At^q — 5a^ . Also: Ar.ah = and [Ar;^] = Ar^a 

[Ar;a];h = AT;afc- 

Turning to the computation of two covariant derivatives: 



Ar., 



= Ar;a,fc - T^ab^T-p 

= Ar;a,fc - r^'ah Ar.p 

— ^ ab ^p 



and three covariant derivatives: 



AT;a 



fcc 



-L a 



fc;c 



Using these results, the expansion tensors ( D2c ) become 

ST^^\b = 5 J 5. 

ST^^\be - r\tSe^ 
5T^^\bcd 



ab ^ cd ^ ab\e^d 



(D3) 
(D4) 

(D5a) 
(D5b) 

(D5c) 
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APPENDIX E: SERIES EXPANSION OF HADAMARD FORM 

For the regularization of the coincident limit of the noise kernel, we need the Hadamard form 

S{x, y) = J^;^ ( — 1- {vo + crt^l + o-^W2) logcr + {awi + (T^W2) ) 



(El) 



to fourth order in cr". We review the standard techniques for finding these expansions and present the results we use. 
The functions Vn{x,y) and Wn{x,y), n > 1 are determined by demanding 

{n~ R/6)S{x,y) ^0. (E2) 

The arbitrary function wo{x,y) is assumed to vanish. Working to fourth order, we proceed in the now familiar 
pattern of expanding each of the function in a s eries expansion and then solve for the expansion tensors by putting 



the expansion in the equations derived from (E2). Using the differential operators 



6 

we need to solve 



Hn = a-PS/p + (^(n - 1 + i (□d))^ (E3a) 
K = D-^ (E3b) 



voix, y) - vT + ^y'v^o^ + ^^^'^Si + '^W'-^g) , + aVVV^T;^?,,, (E4a) 
HqVq = -Kl^l'^ (E4b) 



V, {x,y) « v^f^ + aPv['J + a^a^v^X (E5a) 
2Hivi = -KvQ (E5b) 



V2{x,y) w ^2°^ (E6a) 
m2V2 = -Kvl (E6b) 



along with 



wiix, y) « u.[°^ + cjPw^^^ + aPfj'zwgl, (E7a) 

2HiWi + 2H2vl = (E7b) 

W2{x,y) w (E8a) 

4H2W2 + 2H4V2 = -Kwl (ESb) 

Proceeding as in Appendix we find 

= (E9a) 

= (E9b) 

= {R.ab - 3 Rab-i/ + 4 Rpa Rb" - 2 Rpq Ra'b" + 2 Rpqra Rb"") /360 (E9c) 



'^0 



,(1) 



(2) 
"Oab 



^Oabc ~ R;abc ~ Rpa^bc" ~ Rab^pcF ^ Rab;p^ c ~\~ ^ Rab;cp 
+ 10 Rpa-b R/ + 15 Rab;p RJ' — 10 Rpaqb;c R^'' 
24 Rpq\a Rb^ ~t" 4 Rpa;q Rb^ 8 Rpaqr\b Rc^^ 

-\-2R J, /? P"^*" - 2 7? J, /? + 2 /? I, R ''P'^ 



P 
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+2iip„q^;6ii/^«)/1440 (E9d) 



^otL = {.-^l^GRpa-.bRcd'" -IIQbRab.pRcd'" -™Rpa,bR/;d 

+ 11160 i?pg;a i?fc^«;d - 3600 + 270 i?paq6;r ^ 

-|-2556 Rpqra;b Rc^ d H" 36 Rpqra\b Rc^^ ;d 522 Rpaqb;r 
-522 i?p„,6;^ iJc^^^d + 468 Rpgra;b R/"" ' 1170 ii;„6cd 
+ 1656 i?„a 

•.bed + 1404i?p„;b/<i + 828i?a6 

+702 Rab 810 Rab:p''cd + 828 Rab: cpd 

— 1188 Rab-.cp^d — 1440 Rab-.cdp^ + 210 R;ab Red 

-630 i?a6;pP Red - 3996 i?,p„;5, i?,/ - 4914 Rab.pc Rd" 

+ 1008 Rab;cp R/' + 840 i?pa Rbc R/ + 2916 i?pagfc;cd R"" 
—336 i?;pa Rbcd^ + 1008 Rpa-q' Rbcd^ + 1344 _Rpg _Ra'' Rbcd^ 
—336 i?,pa Rqrb^ Rcd'^^ + 10152 Rpq-^ab RJ' d'^ + 2088 Rpa;qb Rc'd'^ 
-4320 i?pa:hg iJ/d" - 1728 i?a6;pg i?// + 672 i?,p, Rqt R/d" 

-420 i?p, i?// + 4728 i?,p„ i?,qfc/ R,\f + 1536 i?p, iJ.ab'' i?c«/ 

+64 Rpqra Rsb^' Rc^ d'' + 64 Rpqra Rs^ b -Rc^d" + 3168 Rpqra Rs^b^ Rc'd^ 
\.'^^2 Rpqj-Q Rg}^^ Rq fi -\- 72 Rpqj-Q^-i)f^ Rfi^*^ 540 /^pagb;rc -^d^^ 

— 1296 Rpaqb-cr Rd^'^^ — 672 i?pq Rrabc Rd^''^ + 1344 Rpqra Rsbc^ Rd^^^ 
-540 i?paqh;rc i?d«P'' " 1296 Rpaqb:cr Rd'^'"' - 252 Rpa Rqbrc Rd"'"' 
+432 Rpqra-bc Rd^^'^ + 420 i?a6 Rpqrc Rd^^'^ — 252 i?pa Rqbrc Rd 

+ 1344 Rpqra R.sbc^ Rd^^^ — 128 Rpqra Rabc' Rd^'^^ + 96 Rpqra RsbJ Rd'^^'^ 

— 128 Rpqra Rsbc' Rd"'^^ — 672 Rpqra Rsbcd R^^^^ — 1200 Rpaqb Rrcad i?*'^^* 

-624 Rpaqb Rrcsd R^'n /907200 (E9e) 



v^^> = {R.pP - RpqRP^ + RpqrsR^'"'') /360 (ElOa) 

v[^J = {-6 R,paP-UR;pP a + 28 Rpa.q^P- 22 R,pRaP 

-|-40 Rpq-a R^^ 56 Rpa;q R^^ 37 Rpq-r Rc^^ 19 Rpq;r Rq!^^ 

-17Rpqrs;a i?^"" + 2 i?p<,^«;a RP'""- + 12 i?p„g^;, i?f /4320 (ElOb) 

= (210 i?;pi?a6'P + 380 i?;pi?„P;6 - 360 i?p„;,i?5f'« - 2280 i?pa;gi?fc«'^ 

-800 i?p, ; a i?^^ 6 - 200 Rpa:qRP'':b - U88 R.pab'' + 528 i?;p/6 
+820 R;pPab + 180 R.abp" + 960 i?pg;abP« + 960 Rpa^qb""" 
+960i?pa;q6^^ - 1056i?pa;/6f - 944 Ep^jg^Pj + 960 i?p„;f,/9 
-1560 i?pa:fc,«f +480 i?,,,,,;,,/" - 780 Rab-W+ ^80 Rab:p%' 
+ 140 i?;pf i?afe - 152 R.,pa Rb" - 264 i?pa;g« - 1232 i?p,;afa 

+1768i?p„;5bi?f3 - 1056i?p„;b,i?f« - U2Rab;pqR'"' 

-3680 i?pa i?f « - 140 Rpq Rab R""^ 

+1320 Rpq.,r Ra^'b"''^ + 360 Rpq,r Ra^b'^'" + 5560 i?^,;, 

p p pqr-.s i ion p p prs;q i lo p p rpq\s 

(JIU ^pqra;s ^b ' -LOU J^pqra:s ^b i ^pqra;s ^b 

+ 180 Rpqra;s Rb^^^"'^ + 705 i?,p,rs:a R^'^'^^-b " ^716 Rpqra;s R^'^^^;b 
-1008 i?p,,a;s + 544 i?,;p, - 72 i?pq^/ i?,^ fc« 

+1256 Rpq RrP Ra'^b + 6832 Rpq.,ra Rb""'' - 672 i?p,;,, RyP'^'' 

+ 1728 iipa;9r " 4712 Rpq Rra Rb'"''' + 720 Rpa;qr Rb"'"' 

—246 Rpq Rrsa' Rb'^^^ ~ 784 Rpq Rraa^ Rb'^^^ — 276 Rpq RrsaF Rb^'^^ 

-464 i?p, Rras" Rb"' + 1472 i?p,i?,„/i?6-''"^ + 664 i?p,,,;,bi?f«" 

-576 i?p,ra;.6 - 1416 Rpqra-M R'""' ' + 140 Rab Rpqra R'""'' 

— 1020 Rpqra Rstb^ JlPI'^'t ^ 54 ij^^^^.^^ _RP' — 96 Rpqra;ab R^^^^ 
—432 Rpqra:bs R^^^^ — 192 Rpaqb;rs R^^^^ — 2576 i?,,,, Rrasb i?^'"'^* 
—432 Rpqra Rstb'' R^^ '^^ + 1624 _Rpa Rqrsb -R^"'' + 240 Rpqra Ratb^ J^P'^l't 
+240 Rpqra Rstb"^ J{Psrt ^ gg2 ^^^^^ -Rtab'' -R^*' — 16 Rpqra Ratb'' -R*''** 
—2068 Rpqra Rsbt^ i?'"*^ — 1972 Rpqra Rstb^ R'^ 

-2988 i?p,^a i?sfet'' R"*"^') /604800 (ElOc) 



<^ = (I39i?;pi?'f - 354i?pq;^i?f«'''+732i?pg;^i?P'^;« + 114i?;pfq« 
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-288 Rpq./'^^ + 508 Rpq - 84 Rpq./ R^" 
+484 RpqRr^ R''' ) /30'2400 

+ (405 Rpgrs;t R""'''''' + 810 Rpgrs;t R""''*''' + 288 Rpq,rs R""""' 

-5520 Rpq Rrs R^''"' - 2352 Rpq Rrs,^ R"'''" 

-3520 Rpqrs Rt^u R"*'"" - 640 Rpqrs Rtu^" i?'''*") /3628800 (Ell) 



w^^' = - {R.pP - Rpq RP" + Rpqrs R^"''') /2i0 (E12a) 

w['J = {UR.paP + 31R.pPa-?>6Rpa-q'"' + UR.pRaP 

-86 Rpq..aRP' + 112 Rpa-qRP" + 74 Rpq.rRa'"''^ + 5SRpq;rRa'''"' 

+^ORpgrs;a i?^^" " 4it!j,g„;„ E^'^*'' - 24Epa9r;« i?^'^"^) /6480 (E12b) 
"^ll = (-630 i?afc'^ - 1756 R,p RaP;b + 1080 Rpa;q RbP'" + 6840 i?pa;9 Rb"''" 

+3856 i?f «;6 - 968 Rp^.q R^^-^b + 4464 i?.p„t,^ - 1752 i?.^/;, 

-3020 „6 - 540 i?.afc/ - 2880 Rpq-.ab'"' - 2880 il'pa;g6^'' 

-2880 i?pa;9h^P + 3168 i?pa;9'6^ + 3616 Rpa-^q'^b " 2880 it-pa;;,/" 

+4680Ep„;6g9^' - 1440i?„f,;p/« + 2340E„f,;p,«f - 1440i?„b;/,« 
-504 R.,pP Rab - 160 R;pa Rb" + 792 i^p^;," iJ^f + 5152 Rpg.ab R"" 
-6872 i?f 9 + 3168 i?^,.,^ + 395 r^^.^^ rpi 

+ 11040 i?pa RqbRP'^ +504:Rpq RabR^"^ 

-3960 Rpq^rRa'^b"'"' - 1080Rpq..rRj'b'"' - 182A8Rpq;rRj'''^;b 
+ 1548 Rpqra:s i?6^«"'' - 540 Rpqra.s Rb'"''"' - 36 Rpqra.s Rb'""''' 
-540 Rpqra.s Rb""''" - 2955 Rpqrs-.a R'''''':b + 5484 Rpqra;s 

+3024 ii-pg.a,, RP'"'.b - 1632 7?;^, i?/^'' + 216 Rpq./ RJb" 
-3768 Rpq RP Ra^b - 22064 Rpq.ra Rb'"'"' + 2016 7?^,.^^ Rb^"'' 
-5184 iifeP"' + 14136 Rpq iJ^^ -^2160 Rpa.qr Rb""'' 

+738 i?pq Rrsa' Rb'^^'^ + 2352 _Rpq Rras^ Rb''^" + 828 Rrsa' Rb^'^^ 

+ 1392 Rpq Rras^ Rb'^'^'^ — 4416 Rpq Rras^ Rb^'^'^ — 2804 Rpqrs;ab RP'^^^ 
+2Q&ARpqra,sb R""''' + A248 Rpqra;bs R""''' ' 504 R^b Rpqrs -R^'" 
+3060 Rpqra Rstb R""'* - 248 Rpqrs;ab R'"'"' + 288 Rpqra;sb R""""' 
+ 1296 Rpqra:bs R'"'"' + 576 Rpaqb:rs R'"'"' + 7728 Rpq Rrasb R""""' 
+ 1296 Rpqra Rstb' - 4872 Rpa Rqrsb RP'"'' - 720 Rpqra Rstb'' R"'"' 

— 720 Rpqra Rstb"^ i?^"*^' — 2976 Rpqrs Rtab"^ RP*^"^ + 48 Rpqra Rstb^ i?*'^'** 
+6204 Rpqra Rsbt" i?"'"* + 5916 Rpqra Rstb" R"'''' 

+8964 Rpqra Rsbt^ i?^*") /1451520 (E12c) 



= - {139 R.pR'P -354 Rpq,rR^'^'''' + 732 Rpq.rR^'"''' 

+114 R,pPq'i - 288 Rpq.J'^P + 508 R,pq RP'^^ 
-84i?p,./ RP'^ + 484 Rpq Rr^ R"') /145152 

- (405 Rpqrs:t + 810 Rpqrs;t R"'"'''' + 288 Rpq;rs R"""" 

— 5520 Rpq Rrs RP^"^^ — 2352 Rpq Rrst^ i?'*'"* — 3520 Rpqrs Rt^u^ Rqtsu 

-640 Rpqrs Rtu"'^ i?'''*") /1741824 (E13) 



[1] N. D. Birrell and P. C. W. Davies, Quantum fields in curved space (Cambridge University Press, Cambridge, England, 
1982). 

[2] S. A. Fulling, Aspects of quantum field theory in curved spacetime (Cambridge University Press, Cambridge, England, 
1989). 

[3] R. M. Wald, Quantum field theory in curved spacetime and black hole thermodynamics (Cambridge University Press, 
Cambridge, England, 1994). 

[4] A. A. Crib, S. G. Mamayev and V. M. Mostepanenko, Vacuum Quantum Effects in Strong Fields (Priedmann Laboratory 
Publishing, St. Petersburg, 1994). 



21 



[5] B. L. Hu, Physica A 158, 399 (1989). 

[6] E. Calzetta and B.L. Hu, Phys. Rev. D 49, 6636 (1994); B. L. Hu and A. Matacz, Phys. Rev. D 51, 1577 (1995); B. L. Hu 
and S. Sinha, Phys. Rev. D 51, 1587 (1995); A. Campos and E. Verdaguer, Phys. Rev. D 53, 1927 (1996); F. C. Lombardo 
and F. D. MazziteUi, Phys. Rev. D 55, 3889 (1997). 

[7] A. Campos and B. L. Hu, Phys. Rev. D 58 (19 98) 125021: Int . J. Theor. Phys. 38 (1999) 1253. 

[8] B. L. Hu, Int. J. Theor. Phys. 38, 2987 (1999); |gr-qc/9902064 . 

[9] R. Martin and E. Verdaguer, Phys. Rev. D 60, 084008 (1999). D 61, 124024 (2000). 
[10] B. L. Hu and E. Verdaguer, "Recent Advances in Stochastic Gravity: Theory and Issues". Invited Lectures at the Erice 

School, May 2001, in Adv ances in the In terplay between Quantum and Gravity Physics edited by P. Bergmann and V. De 

Sabbata, (Kluwer, 2002). |gr-qc/0110092| 
[11] Ya. Zel'dovich and A. Starobinsky, Zh. Eksp. Teor. Fiz 61, 2161 (1971) [Sov. Phys.- JETP 34, 1159 (1971)]; L. Grishchuk, 

Ann. N. Y. Acad. Sci. 302, 439 (1976); B. L. Hu and L. Parker, Phys. Lett. A 63, 217 (1977); B. L. Hu and L. Parker, 

Phys. Rev. D 17, 933 (1978); F. V. Fischetti, J. B. Hartle and B. L. Hu, Phys. Rev. D 20, 1757 (1979); J. B. Hartle and B. 

L. Hu, Phys. Rev. D 20, 1772 (1979); ibid. 21, 2756 (1980); J. B. Hartle, Phys. Rev. D 23, 2121 (1981); P. A. Anderson, 

Phys. Rev. D 28, 271 (1983); ibid. 29, 615 (1984). 
[12] D. J. Gross, M. J. Perry and L. G. Yafle, Phys. Rev. D 25, 330 (1982). 
[13] B. F. Whiting and J. W. York, Jr. Phys. Rev. Lett. 61, 1336 (1988) 
[14] N. G. Phillips and B. L. Hu, Phys. Rev. D 63, 104001 (2001). 
[15] C.T. Kuo and L. H. Ford, Phys. Rev. D 47, 4510 (1993). 

[16] Chun-Hsien Wu, L. H. Ford, "Quantum Fluctuations of Radiation Pressure" Phys. Rev. D64 (2001) 045010 

[17] N. G. Phillips and B. L. Hu, Phys. Rev. D 62, 084017 (2000). 

[18] N. G. Phillips and B. L. Hu, Phys. Rev. D 55, 6132 (1997). 

[19] Guido Cognola, Emilio Elizalde, Sergio Zerbini Phys. Rev. D 65 (2002) 085031 



[20] Klaus Kirsten Spectral functions in mathematics and physics hep-th/0007251 



[21] R. Camporesi, "Harmonic Analysis and Propagators on Homogeneous Spaces" Phy. Rep 196 (1990) 1-134 
[22] B. S. DeWitt, Dynamical Theory of Groups and Fields, Gordon and Breach, New York, 1965). 
[23] S. M. Christensen, Phys. Rev. D 14, 2490 (1976); ibid. 17, 946 (1978). 

[24] S. M. Christensen, Invited Talk at the International Workshop on Heat Kernel Techniques and Quantum Gravity, University 
of Winnipeg, Canada, August 1994. Proceedings published as Discourses in the Mathematics and Its Apphcations No. 4, 
edited by S. A. Fulling (Texas A & M University Press, College Station, 1995) 

[25] S. L. Adler, J. Lieberman and Y. J. Ng, Ann. Phys. (N.Y.) 106, 279 (1977). 

[26] R. M. Wald, Phys. Rev. D17, 1477 (1977). 

[27] J. D. Bekenstein and L. Parker, Phys. Rev. D23 2850 (1981). 

[28] D. N. Page, Phys. Rev. D25, 1499 (1982). 

[29] B. L. Hu and N. G. Phillips, Int. J. Theor. Phys. 39, 1817 (2000); ^r-qc/0004006[ 



[30] L.H. Ford, Chun-Hsien Wu, "Stress Tensor Fluctuations and Passive Quantum Gravity Int. J. Theor. Phys. gr-qc/0102063 
[31] J. Schwinger, Phys. Rev. 82, 664 (1951); B. S. DeWitt, Phys. Rep. 19C, 297 (1975) 
[32] P. R. Anderson, W. A. Hiscock and D. J. Loranz, Phys. Rev. Lett. 74, 4365 (1995). 
[33] W. A. Hiscock, S. L. Larson and P. R. Anderson, Phys. Rev. D (1997) 

[34] L. Parker and S. M. Christensen, MathTensor: a system for doing tensor analysis by computer (Addison- Wesley Publishing 

Company, Redwood City, California, 1994). 
[35] S. Wolfram, Mathematica: a system for doing mathematics by computer (Addison- Wesley Publishing Company, Redwood 

City, Cahfornia, 1991). 

[36] We must stress, as we did in Paper I, that though all these are done on a computer, no numerical approximation is used. 

All work is done symbolically in terms of the explicit functional form of the metric and the parameters of the field. The 

final results are exact to the extent that the analytic form of the Green function is exact. 
[37] Since this work for analyzing the coincident limit of the noise kernel via the point separation method is tailored to symbolic 

computation on the computer, our method for deriving particular results is designed to take maximum advantage of the 

computer. 



